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ABSTRACT 

Precision surface finishes are used in a wide variety of applications. From bearing races and rolling elements 
to parallel slide ways, the topography of these surfaces is critical to the performance of the products. With the 
advent of diamond turning and other abrasive precision finishing such as magneto -rheological finishing 
process, nanometric surface finish can be realized. The ultra precision surfaces also find applications in optics 
and/or bio-medical implant surfaces. The 2D roughness parameter defined in standards are inadequate for 
characterising surfaces therefore analysis of 3D roughness parameter is necessary for broad analysis of surface 
conditions. 
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I. INTRODUCTION 

The interaction between two surfaces is an important area of study in industrial applications. A lot of 
research has been done on the surface topography as it helps in understanding the interaction between the 
surfaces. The two-dimensional (2D) parameters defined in standards are inadequate for characterizing surfaces 
because real surfaces interact in 3 dimensions. It can be seen in case of two surfaces in contact 3-D amplitude 
and spatial parameters are important. Also, texture has important role in hemodynamics (dynamics of blood 
flow in and around a cardiovascular implant. The project aims to characterize select engineered precision 
finished surfaces and/or textured surfaces to improve functional response. 




Figure 1.1 The surface interaction in contact of two surfaces. 



Therefore a set of 3D parameters are required to describe the properties of surface in effective way. 
Also the functional properties cannot be obtained by using 2D parameters. Various 3D parameters have been 
defined based on different properties like amplitude properties and spatial properties required for describing the 
surface [1]. These parameters have been successful in analyzing various problems in industrial and research 
field. 

II. 3-D SURFACE PARAMETERS 

2.1 3-D Surface parameters. 

Until recently only 2-D roughness parameters [1-2] were in use which do not provide the complete picture, 
especially concerning the lay characteristics of the surface. Hence, there has been a push within both industry 
and academia for transition to 3-D parameters. This need arises partially due to the realization that surfaces 
interact in three dimensions instead of two [3]. The primary parameter set proposed by Stout et al. [7] involves 
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14 parameters, which characterize the major aspects of topographic features. There are four parameters for 
describing the amplitude and height distribution properties, four parameters for describing spatial properties, 
three parameters for describing hybrid (i.e. both amplitude and spatial) properties and three parameters for some 
functional properties. Table 2.1 shows the 14 3-D parameters proposed by Stout et al [7]. 



Table 2.1 - Set of 14 3-D parameters given by Stout et al. [4] 



Amplitude Parameters 


Spatial 
Parameters 


Hybrid Parameters 


Functional 
Parameters 


RMS Deviation 

s q 


Density of 
Summits 

Sds 


RMS Slope 
SAq 


Surface Bearing Index 

Sbi 


Ten point Height 

S z 


Texture Aspect 
Ratio, 

Str 


Mean Summit 
Curvature 

S S c 


Core Fluid Retention 
Index 

S c i 


Skewness 

S s k 


Texture 
Direction 

Std 


Developed Area Ratio 

Sdr 


Valley Fluid Retention 
Index 

S v i 


Kurtosis 

Sku 


Fastest Decay 
Autocorrelation 
Length 

Sal 







2.2 AACF and APSD Analyses 

The autocorrelation function for surface profiles is based on the idea that profile readings can be treated 
as random signals. This type of analysis has found extensive use in random signal processing in communication. 
The autocorrelation function for surface profiles is based on the idea that profile readings can be treated as 
random signals. This type of analysis has found extensive use in random signal processing in communication. 
The idea of computing various random process functions for machined surfaces was first proposed by Dong et 
al. [7] used AACF and APSD for characterization of 3-D spatial parameters. They also studied characteristics of 
AACF and APSD for a variety of engineered surfaces. The following section provides a detailed description of 
AACF and APSD which form the basis of 3-D spatial surface characterization. 

2.2.1 AACF 

The areal autocorrelation function describes the general dependence of data at one position on the data 
at another position. As discussed earlier, the autocorrelation function has been effectively used for random 
signal processing. It provides information about time data in case of communication signal and spatial data in 
case of 2-D surface analysis. For 3-D surface analysis the areal ACF or AACF is given as follows: 

R(t x ,t y ) = E[tj(x, y)rj(x + r x , y + r y )] 

l y I 
1 X 

= lim f f 77 (jc, y )tj (x + r ,y + r )dxd y (2.1) 

4/ / J J ^ 

x y ~ l y -K 

A non-biased discrete estimation of the AACF is given by Eq. (2.2). 

* N - j M - i 

R(*i>*j) = Z Z ^K^z^K^') 7 /^) (2.2) 

(M -i)(N - j) l=l k __ x 

where, 

i = 0,1...., m < M ; j = 0,1...., n < N 
t . = /Ax; t = jAy 

rj (x, y) is the residual surface [5] after a plane is fit to remove the form (longer wavelength and undulations) 

from the surface data; m and n are the autocorrelation lengths in the x and y directions, respectively. As i, j 
approach M and N, respectively, fewer data points are available for computation of AACF and the statistical 
confidence decreases considerably. Therefore, maximum coordinates in the x and y axes are restricted to 
x m =AxM/2 and x n =AyN/2. Figure 2.1 shows the normalized AACF plots for the surfaces given by Singh et al. [5] 
for ground (GD), honed (HN), hard turned (HT) and isotropic finished (IF) surfaces used in bearing industry. 
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Figure 2.1. Normalized AACF plots for the: (a) GD, (b) HN, (c) HT and (d) IF surfaces [5]. 



The AACF is an effective tool to visualize the correlation in different directions. It is evident from 
Figure 2.1 (a) that for the GD surface, which exhibits a strong directional property, the decay of the auto 
correlation function is very slow along the lay and very rapid across the lay. For HT surface (Figure 2.1(c)) a 
distinct periodicity can be observed. For the IF surface, there is no correlation in any particular direction and 
hence the decay is uniform in all directions as shown in Figure 2.1(d). Thus, it can be inferred that distinctive 
properties can be identified from the AACF. Consequently, AACF analysis is used to compute 3-D texture 
parameters, namely texture aspect ratio (S tr ) and fastest decay (S a i). 

2.2.2 APSD 

Spectral analysis is a powerful tool to condense the time/space domain information into the frequency 
domain. Once frequency domain data is available we can accentuate the information for individual frequency or 
wavelength contribution. Power spectral analysis can help in distinguishing between the roughness and 
waviness because of distinct differences in the wavelengths. Dong et al. [8] discussed the Hermitian symmetry 
for 2-D spectral analysis and algorithm to compute the APSD and extract texture parameters. For the residual 
surface r|(x,y) (given by Dong et al. [9]), its continuous and discrete Fourier transforms are given by Eqs.(2,3) 
and (2.4), respectively. 

00 00 

F (m w )= \7](x,y)e dxd y (2.3) 



N-1M-1 ' 2 /T ( k \ / ) 

^ P ,*%) = I 2><** + i.:v, + i)« " N (2-4) 

1 = 1 k = l 

where m and w are angular frequencies in two directions as defined below: 
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p = 0,1 M - 1; q = 0,1 N - 1; 

G7 = , tU — 

P AxM q AyN 
The areal PSD can be implemented by: 

G (m , m ) = F (m ,m )F (m , m ) (2.5) 

K p 7 q ' ^p 1 q'^p 1 q / v/ 

where '*' denotes complex conjugate. Manipulations need to be carried out to get the center coordinates to be at 
the origin. Figure 2.2 shows the APSD plots for the various surfaces [5]. 

For anisotropic surfaces, most of the power is located in the direction perpendicular to the lay. This is 
observed in the GD surface (Figure 2.2 (a)) and to some extent in the HN surface (Figure 2.2 (b)). The HT 
(Figure 2.2 (c)) surface has power concentrated at the specific frequency of periodicity perpendicular to the lay. 
The IF surface (Figure 2.2 (d)) is a random surface and hence the power is concentrated at the origin. The APSD 
is required to compute the texture direction S td parameter. 




q - 1 

(c) (d) 

Figure 2.2. APSD plots for the: (a) GD, (b) HN, (c) HT and (d) IF surfaces [5]. 



2.3 Computation of 3-D Parameters 

As discussed previously, there exist a set of 3-D parameters that describe the spatial and amplitude 
features of a given surface. These parameters provide a more realistic estimate of the 3-D surface features in 
comparison to their 2-D counterparts. Some of these 3-D parameters were computed for the four surfaces of 
interest here. These parameters are discussed in detail by Dong et al. [3-4]. 



2.3.1 3-D Amplitude Parameters 

Amplitude parameters are important parameters for surface topography. Many 2-D parameters [1-2] 
have been defined in national and international standards. The two 3-D amplitude parameters that were 
computed for the surfaces studied here are: 3-D root mean square roughness (S q ) and surface skewness (S sk ). 
Root Mean Square Deviation of Surface Topography (S q ) 

This is a statistical amplitude parameter defined as the root mean square value of the surface asperity departures 
from the reference datum within the sampling domain. This is assumed to provide a conservative estimate of the 
average asperity height measured from the datum. The S q parameter, in continuous and discrete forms, is given 
as follows: 



S = 



N M 



1 ff 2 I 1 

jj77 (x, y)dxd y = J— — ^ J] 1 (x.,y.) 



(2.6) 



/ / 

x y 0 0 



MN 



7=1 i=\ 
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where, l x and l y are the side lengths of the sampling area. This parameter is of significant tribological 
importance as it plays a key role in determining the factional response of a surface as will be seen later in the 
report. 

Skewness of topography height distribution (S sk ) 

The skewness is the measure of asymmetry of surface deviations about the mean plane. It is given defined as 
follows: 



N M 

3 



1 3 1 

S = [ \r/ (x, y) p(ri)dxd y = S 77 K'^P ( 2 ' 7 ) 

S MNS /=1 i- x 

q -00-00 q J 1 * 1 

It can be used effectively to describe certain aspects of the shape of the topographic height distribution. From a 
functional point of view, it can indicate presence of spikes. Large negative skewness may be found in surfaces 
that have pits and troughs. S sk is unstable if used for indicating functional properties of the surface because the 
effect of a small number pits or spikes can have a significant effect on the computation of skewness. 

2.3.2 3-D Spatial Parameters 

Spatial parameters are important because 2-D parameters cannot capture directionality and the degree 
of anisotropy. Both deterministic and random surfaces can possess isotropic or anisotropic character. The 
greatest difficulty with the spatial parameters is that they are usually difficult to describe by parameters owing to 
the random and multi- wavelength components of the surface. The four important spatial 3-D parameters 
evaluated here are: S a i, S tr , S td and S ds 
Fastest Decay Autocorrelation Length (S a i) 

This parameter describes the nature of the AACF. It is the shortest distance from the origin in which the 
normalized AACF decays to a value of 0.2. Whitehouse and Archard [10] modeled the autocorrelation decay for 
a random surface with Gaussian distribution to be exponential. The fastest autocorrelation decay to a value 0.1 
was used. Since actual surfaces do not display exponential decay, Dong et al. [4] used a decay value of the 
normalized AACF of 0.2. Mathematically, S a i is given as follows: 



S„, =min( J*) +r 2 y \R{T x ,T y )< 0.2 (2.8) 



This is an important parameter as it identifies the direction in which the correlation is less. It decays very rapidly 
across the lay direction for a surface exhibiting strong lay characteristics as illustrated in Figure 2.1 (a). Also, 
the decay rate of the AACF is governed by the significant frequency components. A large value of S a i indicates 
that the surface is dominated by low frequency components whereas a small value indicates otherwise. Turned 
surfaces have larger values and ground and honed surfaces have smaller values of S a i. 
Density of Summits (S ds ) 

Density of summits is the number of summits contained in a unit of sampling area. It is computed as follows: 

Number of Summits 

S ds = (2.9) 

(M -l)(N - l)AxAy 



This parameter is a function of summit definition. A summit is a central point of an area above the 
mean plane which has highest value among all data points of the area. The definition of sampling area is 
ambiguous; there are many definitions that exist such as the nearest eight neighbors [11], a contour area or a 
window. Thomas [11] showed that the S ds (based on 8 nearest neighbors) increases rapidly as the sampling 
frequency increases. This necessitates an ideal definition of area over which a summit is defined. Dong et al. [4] 
suggested a sampling area based on S a i and S ds computed from this definition is not very sensitive to sampling 
frequency. Sampling area is a square (Ax=Ay) whose half length is equal to the S a i. Although more summits can 
be found using any other definition, they are considered correlated hence only one summit is recognized. In this 
study S ds is computed based on S a i. The abrasive processes such as the ground and honed surfaces have very 
high frequency components, meaning that the lay is very closely spaced and therefore summits occur more 
often. Consequently, the S ds value is higher for those processes. It is a very good estimate of the average number 
of asperities per unit area and is found to influence the frictional response of the surfaces under mixed 
lubrication. 

Texture Aspect Ratio (S tr ) 

The texture aspect ratio is used to identify the texture pattern, i.e. long-crestedness or isotropy. Texture 
aspect ratio can serve as the measure of isotropy. The definition of texture aspect ratio is as follows: 
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s„ = 

The distance that the normalised AACF has the fastest decay to 0.2 in any possible direction 
The distance that the normalised AACF has the slowest decay to 0.2 in any possible direction 

(2.10) 

A large value of this parameter indicates strong isotropy (S tr > 0.5) while a small value (S tr < 0.3) 
indicates strong anisotropy. In the study done by Singh et al. [5-6], IF surface has a value of 0.67. The GD 
surface, which is strongly anisotropic, exhibits a very low value of 0.08. 
Texture Direction (S td ) 

Texture direction is used to indicate the pronounced direction of lay. This is extracted from the APSD. Since the 
coordinate system is arbitrary, the coordinate system [4] is defined as shown in Figure 2.3. The formula for 
computation of S td is given by: 




n 



(2.11) 



= n - p , V < p < n 



where p= value of 0 at which G a (0) is maximum. 




Figure 2.3. Definition of texture direction. 



G a (0) is the angular spectrum it is derived from APSD by integrating the spectral energy radially between 0 and 
179 degrees. 

G a (0)= jG(0,r)dr 



R(0) = -[(Ax cos 6>) 2 + (Ay sin <9) 2 ]" 1/2 
2 

0 < 6 < 179 



(2.12) 



First G(0,r) was extracted from the APSD for all the surface orientations and then sorted in ascending order. 
Numerical integration was carried out by the trapezoidal rule and the maximum value computed. 





o x us as ao iao la 140 la leo xo 



0 20 40 60 60 1C0 120 140 160 180 200 



0 (Degrees) 



0 (Degrees) 



(a) (b) 

Figure 2.4. Typical angular spectrum plot for texture orientation (a) lay perpendicular to Y-axis and (b) isotropic 

surface 
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The texture direction is given by the angle corresponding to the maximum amplitude in the angular 
spectrum plot. Since all surfaces except the IF surface have lay perpendicular to the Y axis, the maximum 
amplitude occurs at an angle of 90° as illustrated in Figure 2.4(a). The IF surface has no preferred direction of 
lay and consequently there is no peak in the angular spectrum plot as shown in Figure 2.4 (b) [5]. 



2.3.3 Hybrid Parameter 

RMS Slope( S Aq ) 

The root mean square slope of the surface within the sampling area is computed as follows: 



r 



V 



tj(x., y .) - 77 (x._ t , y .) 



Ax 



J 



+ 



V 



7] (x., y .) - ri(x._ x , y . ) 



Ax 



\ 



J 



1 



1/2 



(2.13) 



N M 



0.5 



I5> 

7=2 i=2 



2 
U 



Aq 



(M - l)(N - 1) 



(2.14) 



V J 

Arithmetic mean summit curvature of the surface (S sc ) 



III. CONCLUSION 

[I] Experience shows advantages of certain finishes in certain roles. 

[2] 2D roughness measurement(stylus) often correlated with functionality. 

[3] Surfaces contact in three dimensions. 

[4] 3D texture parameters affect functionality. 

[5] AACF describes the general dependence of data at one position on the data another position used for fastest 
decay autocorrelation length and density of summits. 

[6] To calculate the spatial parameters, the AACF, APSD and Angular spectra for each surface must be 
analysed. 
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